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RESPIRATION  AND  HEMOLYSIS 
OF  ULTRAVIOLET 

IRRADIATED  FROG  ERYTHROCYTES1 

Lewis  E.  Harvie  and  S.  P.  Maroney,  Jr. 

University  of  Virginia 

INTRODUCTION 

Ionizing  radiations  have  been  known  for  some  time  to  cause  hemo¬ 
lysis  in  mammalian  erythrocytes.  Early  experiments  indicated  that  X- 
irradiation-induced  hemolysis  of  human  and  beef  erythrocytes  was 
preceded  by  slow  swelling  of  the  cells  (Ting  and  Zirkle,  ’40).  In  the 
same  report,  these  authors  suggested  that  hemolysis  and  swelling  were 
dependent  on  an  increased  permeability  to  sodium  and  potassium  ions. 
Subsequent  analysis  by  Sheppard  and  Beyl  (’51)  confirmed  the  hypo¬ 
thesis  that  the  cell  was  more  permeable  to  sodium  and  potassium  after 
X-irradiation;  they  further  showed  that  these  ions  exchange  quantita¬ 
tively.  Non-ionizing  radiations  have  also  been  shown  to  cause  hemolysis 
(Sonne,  ’29;  Cook,  ’56;  Green,  ’56).  The  hemolytic  effects  of  ultraviolet 
light  as  reported  by  Cook  (’56)  and  Green  (’56),  are,  in  general,  in 
agreement  with  those  produced  by  X-rays. 

It  has  been  demonstrated  that  the  energy  for  maintaining  the 
balance  of  potassium  and  sodium  in  non-nucleate  erythrocytes  is  supplied 
by  glycolysis  (Maizels,  ’51;  Weller  and  Taylor,  ’51;  Solomon,  ’52). 
Green  (’56)  has  shown,  however,  that  ultraviolet  irradiation  of  human 
red  cells  causes  an  acceleration  in  the  loss  of  potassium  and  gain  of 
sodium  ions  without  inhibiting  glycolysis.  Sheppard  and  Stewart  (’52) 
obtained  similar  results  with  human  red  cells  exposed  to  X-rays.  This 
suggests  that,  in  addition  to  metabolic  energy,  other  factors  are  impor¬ 
tant  in  regulating  the  exchange  of  these  cations  in  non-nucleated  red 
cells.  One  obvious  structure  for  such  regulation  is  the  cell  membrane 
and,  indeed,  gross  changes  in  the  mammalian  erythrocyte  membrane 


1  This  paper  was  presented  to  the  Virginia  Academy  of  Science  on  May  8,  1959  at  the 
University  of  Virginia.  An  abstract  has  been  published  in  the  Virginia  Journal  of 
Science,  Sept.  1959,  10(4)  :  263-264. 
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have  been  observed  following  X-irradiation  (Zacek  and  Rosenberg,  ’50). 
It  seems  clear,  then,  that  radiation-induced  hemolysis  in  the  mammalian 
erythrocyte  is  preceded  by  a  change  in  ion  permeability  which  is  not 
dependent  on  an  alteration  in  the  energy  supply  of  the  cell,  i.e., 
glycolysis. 

Most  of  the  investigations  on  erythrocytes  have  been  concerned  with 
mammalian  cells;  comparatively  few  studies  on  nucleated  erythrocytes 
are  available.  Hunter  and  Hunter  (’57)  have  shown  that  the  oxygen 
consumption  of  nucleated  erythrocytes  is  significantly  higher  than  that 
of  non-nucleated  erythrocytes  when  compared  at  the  same  temperature. 
Evidence  for  an  active,  oxidative  metabolism  in  nucleated  red  cells  of 
the  chicken  may  be  found  in  the  work  of  Rubinstein  and  Denstedt  (’53), 
and  Defendi  and  Pearson  (’55).  Ion  transport  in  nucleated  erythrocytes 
(chicken,  grass-snake)  depends  on  this  aerobic  respiration  (Maizels, 
’54).  The  problem,  then,  arises  as  to  the  effect  of  hemolvtic  doses  of 
radiation  on  the  aerobic  metabolism  of  the  nucleated  red  cell. 

Martin  and  her  co-workers  (’56),  in  studying  bullfrog  erythrocytes, 
observed  an  inhibition  in  oxygen  consumption  of  these  cells  following 
X-irradiation„  The  present  work  has  been  undertaken  to  study  the  rela¬ 
tion  between  ultraviolet-induced  hemolysis  and  respiration  following  ul¬ 
traviolet  radiation  in  nucleate  erythrocytes  of  the  leopard  frog,  Rana 
pipiens.  The  results  indicate  that  this  radiation  increases  hemolysis  with¬ 
out  inhibiting  respiration. 


METHODS 

Frogs  used  in  this  investigation  were  kept  in  covered  porcelain 
pans  in  about  an  inch  of  water,  the  temperature  of  which  varied  between 
20°C  and  22°C. 

Blood  from  the  aorta  of  one  or  two  pithed  frogs  was  drawn  into  a 
fine  glass  pipette,  the  tip  of  which  had  an  inside  diameter  of  0.7  mm, 
and  mixed  with  phosphate  buffered  Ringer’s  (0.11M  NaCl,  0.02M  KC1, 
0.01M  CaCl2,  0.0009M  Na2HP04,  and  NaH2P04;  pH  7.4)  in  which 
a  few  crystals  of  sodium  citrate  had  been  dissolved  as  an  anticoagulant. 
Just  before  drawing  up  the  blood,  the  pipette  was  rinsed  with  some  of 
this  solution  in  order  to  prevent  clotting  of  the  blood  on  contact  with 
the  inside  wall  of  the  pipette.  Erythrocytes  harvested  from  this  blood  by 
centrifugation  were  washed  twice  and  finally  suspended  in  phosphate 
buffered  Ringer’s  solution.  Cell  concentration,  determined  from  hemacyto¬ 
meter  counts,  ranged  from  70,000  to  150,000  per  cubic  mm. 

The  cell  suspension  was  divided  into  two  3.5  ml  aliquots,  each  of 
which  was  placed  in  a  50  ml  round  bottom  quartz  flask  and  rotated  at 
five  rpm  to  keep  the  cells  in  suspension.  One  aliquot  was  irradiated 
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with  a  general  Electric  G15T8  germicidal  lamp,  the  times  of  irradiation 
ranging  from  10  minutes  to  60  minutes.  The  distance  of  the  surface  of 
the  cell  suspension  from  the  lamp  was  4.5  cm.  All  experiments  were 
carried  out  at  room  temperature  (approximately  25°C).  The  second 
aliquot,  which  was  treated  in  the  same  way  as  the  first  except  that  it 
received  no  ultraviolet  radiation,  acted  as  the  control.  Respiration  of 
both  irradiated  and  non-irradiated  cells  was  measured  using  a  standard 
Warburg  technique,  following  the  procedures  described  by  Umbreit 
et  al  (’57).  Flasks  containing  two  ml  of  cell  suspension,  to  which  was 
added  glucose  in  a  final  concentration  of  50  mg  %,  were  placed  in  a 
water  bath  at  27.3  °C  and  oscillated  100  times  per  minute  through  a  five 
cm  stroke.  Since  cell  number  varied  with  each  experiment,  oxygen 
consumption  was  recalculated  on  the  basis  of  100,000  cells. 

Separate  flasks  containing  about  one  ml  of  cells  were  oscillated  in 
the  water  bath.  For  hemolysis  determination,  samples  were  taken  from 
these  flasks,  diluted  1:10  with  Ringer’s,  and  placed  on  glass  slides  for 
counts  of  hemolysed  cells.  Percent  hemolysis  was  calculated  as  the 
number  of  hemolysed  cells  per  total  cells  counted  (total  cell  counts 
between  400  and  500). 


RESULTS 

Non-hemolysed  cells  had  an  opaque,  yellowish  appearance,  with 
the  nucleus  clearly  discemable.  Dark  bands  running  at  right  angles  to 
the  long  axis  of  the  cell  were  apparent.  These  have  been  described  by 
Trotter  (’56)  who  attributed  them  to  the  residual  alkali  in  the  glass 
slide.  As  the  cell  swelled  following  irradiation,  the  bands  disappeared, 
the  cell  became  subspherical  and  at  hemolysis  rapidly  faded  with  the  loss 
of  hemoglobin  to  a  “ghost.”  The  appearance  of  the  “ghost”  was  the  end 
point  in  these  experiments. 

Cumulative  oxygen  consumption  increased  as  a  direct  linear  function 
of  time  for  about  the  first  six  hours  from  the  beginning  of  the  mano- 
metric  measurements.  A  typical  experiment  (30  minutes  ultraviolet 
radiation)  is  presented  in  Figure  1.  It  should  be  noted  that  there  was 
no  significant  difference  in  the  oxygen  consumption  of  irradiated  and 
non-irradiated  cells. 

Hemolysis  of  irradiated  erythrocytes  increased  in  a  similar  manner 
for  about  the  first  four  hours.  This  was  followed  by  a  leveling  of  the 
hemolysis  curve.  This  “plateau”  in  the  curve  lasted  three  or  four  hours 
after  which  it  rose  sharply  to  100%.  Time  required  for  complete  hemo¬ 
lysis  was  approximately  twelve  hours.  Non-irradiated  cells  showed  no 
hemolysis. 

The  initial  amount  of  hemolysis  after  irradiation  varied  according 
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MICROLITERS  OF  OXYGEN  /  I05  CELLS 


FIGURE  1 

Oxygen  consumption  and  percent  hemolysis  of  ultraviolet  irradiated  and 
non-irradiated  frog  red  cells  (open  circles  =  respiration,  ultraviolet; 
closed  circles  =  respiration,  control;  open  triangles  =  hemolysis,  ultra¬ 
violet;  closed  triangles  =  hemolysis,  control.)  Ultraviolet  dose  =  30 
minutes. 
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MICROLITERS  OF  OXYGEN  /  1 05  CELLS 


FIGURE  2 

Respiration  rates  and  per  cent  hemolysis  values  of  ultraviolet  irradiated 
and  non-irradiated  frog  red  cells  four  hours  after  exposure  to  radiation 
(open  circles  =  respiration,  ultraviolet;  closed  circles  =  respiration, 
control;  triangles  =  hemolysis).  The  control  respiration  rate  is  included 
for  each  ultraviolet  dose  to  better  point  up  the  fact  that  there  is  no 
significant  difference  in  respiration  rate  between  irradiated  and  non- 
irradiated  cells. 
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to  the  dose  used.  At  a  10  minute  dose,  no  hemolysis  was  observed  in 
the  irradiated  cells  until  three  or  four  hours  following  irradiation.  At  a 
15  minute  dose,  hemolysis  began  about  one  hour  after  irradiation, 
while  for  a  20  minute  dose,  hemolysis  began  immediately  following  irra¬ 
diation.  At  30  and  60  minute  doses,  4  and  10  percent,  respectively,  of 
the  cells  were  hemolysed  in  the  first  count  following  irradiation. 

Respiration  rate  for  the  first  four  hours  and  the  per  cent  hemolysis 
at  four  hours'  following  irradiation  are  presented  as  functions  of  ultra¬ 
violet  dose  in  Figure  2.  The  per  cent  hemolysis  increased  with  the 
ultraviolet  dose  while  the  respiration  rate  was  not  significantly  different 
between  the  irridiated  and  non -irradiated  cells. 

DISCUSSION 

The  first  apparent  change  following  ultraviolet  radiation  was  a 
swelling  of  the  frog  red  cell.  A  reasonable  explanation  for  this  swelling 
involves  the  phenomenon  of  colloid  osmotic  pressure  which  consists  es¬ 
sentially  of  two  phases  (Ponder,  ’48):  (1)  a  lytic  agent  alters  the 

cell  membrane’s  selective  permeability  to  cations;  (2)  as  potassium  and 
sodium  ions  diffuse  to  equilibrium  across  the  membrane,  water  enters 
the  cell  under  the  influence  of  the  colloid  osmotic  pressure  of  the  hemo¬ 
globin.  This  penetration  of  water  results  in  swelling  of  the  cell  to  a 
critical  volume  at  which  hemoglobin  is  released  to  the  environment.  The 
phenomenon  of  colloid  osmotic  pressure  has  been  suggested  to  account 
for  hemolysis  following  both  ultraviolet  (Leu  et  al,  ’42;  Cook,  ’56)  and 
X-radiation  (Ting  and  Zirkle,  ’40)  in  mammalian  red  cells.  The  inhibi¬ 
tion  of  ultraviolet-induced  hemolysis  of  frog  red  cells  by  sucrose 
(Maronev,  ’59)  lends  support  to  the  colloid  osmotic  hemolysis  theory. 

For  the  colloid  osmotic  pressure  of  hemoglobin  to  be  effective  in 
causing  hemolysis,  the  ions  must  be  free  to  diffuse  to  equilibrium  across 
the  cell  membrane.  Energy  is  used  to  prevent  the  potassium  and  sodium 
ions  from  reaching  diffusion  equilibrium  in  the  intact  cell  (Harris,  ’41; 
Danowski,  ’41;  Maizels,  ’54).  Inhibition  of  metabolism  is,  then,  one 
means  for  allowing  ions  to  reach  diffusion  equilibrium.  An  example  of 
this  may  be  found  in  mammalian  red  cells  kept  at  low  temperatures 
(5°C).  Potassium  is  lost  from  these  cells  and  regained  upon  returning 
the  cells  to  physiological  temperature,  provided  that  glucose  is  present 
(Harris,  ’41).  The  results  of  the  present  study  show  that  ultraviolet- 
irradiated  frog  erythrocytes  respire  at  the  same  rate  as  non-irradiated 
erythrocytes;  yet  cells  exposed  to  ultraviolet  radiation  exhibited  hemo¬ 
lysis  while  non-irradiated  cells  did  not.  It  seems  apparent  that  the  cell’s 
aerobic  respiration  is  intact  and  an  explanation  other  than  the  inhibi¬ 
tion  of  energy  metabolism  must  be  sought  to  explain  hemolysis  in  the 
frog  red  cell. 
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The  concept  that  metabolic  energy  is  solely  responsible  for  main¬ 
taining  cation  gradients  has  been  criticized  by  Parpart  and  Green  (’51). 
Rabbit  red  cells,  when  treated  with  n-butyl  alcohol,  gain  sodium  and 
lose  potassium  at  comparable  rates.  Further,  the  recovery  of  potassium 
and  sodium  exchange  after  removal  of  the  alcohol  was  much  more  rapid 
than  would  be  expected  if  metabolism  had  been  interfered  with.  They 
point  out  that  it  is  more  likely  that  the  butyl  alcohol  has  an  effect  on 
the  membrane  alone,  permitting  a  rapid  exchange  of  potassium  and 
sodium. 

It  is  suggested,  then,  that  ultraviolet  radiation,  in  the  case  of  the 
frog  red  cell,  altered  the  cell  membrane  allowing  potassium  and 
sodium  to  diffuse  freely.  The  cell  then  swelled  because  of  the  excess 
intracellular  osmotic  pressure.  In  this  regard,  frog  erythrocytes  were 
acting  as  did  the  mammalian  red  cell  following  ultraviolet  radiation 
(Green,  ’56).  It  should  be  noted  that  the  respiratory  rates  of  irradiated 
and  non-irradiated  cells  were  the  same  in  spite  of  the  fact  that  up  to  40 
per  cent  of  the  irradiated  cells  had  hemolysed.  This  would  indicate  that 
the  respiratory  apparatus  of  the  hemolysed  cell  was  intact  and  func¬ 
tioning  normally  although  further  study  of  this  phenomenon  is  desirable. 

SUMMARY 

Frog  erythrocytes  washed  in  phosphate  buffered  Ringer’s  solution 
were  exposed  to  ultraviolet  light  for  various  lengths  of  time.  Respiration 
and  hemolysis  of  both  irradiated  and  non-irradiated  cells  were  measured 
for  about  six  hours  following  exposure  to  ultraviolet  radiation. 

It  was  found  that  this  radiation  increases  hemolysis  without  affect¬ 
ing  oxygen  consumption  of  these  cells. 

The  similarity  of  the  mechanism  of  hemolysis  in  frog  erythrocytes 
to  the  proposed  mechanism  of  hemolysis  in  mammalian  erythrocytes  is 
discussed. 
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BALLOON-WING,  A  NEW  MUTATION 
IN  BLATTELLA  GERMANIC  A  (L.) 

Donald  G.  Cochran  and  Mary  H.  Ross 
Virginia  Agriculture  Experiment  Station 

INTRODUCTION 

The  study  of  spontaneous  and  induced  mutations  in  such  insects 
as  Drosophila,  Habrobracon,  and  several  others,  has  been  of  great  signi¬ 
ficance  in  the  development  of  formal  genetics.  In  fact,  were  it  not  for 
Drosophila,  genetics  might  have  progressed  at  a  much  slower  rate.  With 
regard  to  inheritance  in  the  myriad  of  other  insect  species,  however, 
there  is  still  very  little  information  available.  For  example,  the  mutation 
reported  in  the  present  paper  is,  to  the  writers’  knowledge,  the  first 
genetically  controlled  modification  of  the  wing  to  be  described  in  the 
ubiquitous  German  cockroach,  Blattella  germanica  (L.).  Although  it 
might  be  argued  that  this  is  an  extremely  stable  species  and  that  muta¬ 
tions  would,  therefore,  be  uncommon,  the  facts  do  not  seem  to  bear  this 
out.  During  the  course  of  the  present  investigation  several  other  spon¬ 
taneously  occurring  variants  have  been  found  in  the  same  inbred  stock. 
Some  of  these  traits  are  now  under  study,  and  the  results  will  be 
published  elsewhere. 

In  Drosophila  it  is  known  that  approximately  150  genes  influence 
or  modify  the  development  of  the  wing  (Emerson,  1955:  648).  Muta¬ 
tions  at  several  of  these  loci  show  some  similarity  to  the  one  described 
here.  In  their  compilation  of  known  mutations  in  D.  melanogaster, 
Bridges  and  Brehme  (1944:  20,  23,  29,  26,  152,  95)  listed  the  following 
traits:  balloon  (ba),  Bubble  (Bb),  blister  (bs),  bloated  (bio),  puff 
(puf),  inflated  (if)  and  vesiculated  (vs).  Mutations  of  this  type  are 
also  known  in  other  species  of  Drosophila  (Gordon  et  al.,  1939:  51; 
Metz  et  al.,  1923:  16)  and  in  Musca  domestica  (Sullivan  and  Hiroyoshi, 
1960:  214;  Milani,  1956:  225).  These  mutations  are  all  characterized  by 
swellings  on  the  wing,  usually  filled  with  fluid. 

The  name  “balloon”  was  chosen  for  the  present  mutation  because 
the  character  seems  to  resemble  most  closely  balloon-wing  as  described 
in  the  cheese  skipper,  Piophila  casei,  (Painter,  1917:  306)  and,  to  a 
lesser  degree,  ballon- wing  in  D.  melanogaster.  Balloon  in  P.  casei  and 
in  D.  melanogaster  is  recessive  and  is  not  sex-linked.  As  will  be  shown, 
this  is  also  true  in  B.  germanica. 
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DESCRIPTION  OF  CHARACTER 

Balloon-wing  roaches  are  characterized  by  a  single,  bubble-like 
swelling  on  each  forewing  and,  usually,  in  each  hindwing.  During  the 
course  of  this  study  991  balloon-wing  specimens  were  examined.  On  the 
forewings  the  swellings  are  oval  or  pear-shaped,  usually  from  0.24-0.28 
cm.  in  length  and  0.15-0.17  cm.  in  width.  In  85%  of  the  specimens  ex¬ 
amined  the  anterior  ends  of  the  swellings  occurred  opposite  or  proximal 
to  the  terminus  of  the  first  anal  vein  (1  A),  but  not  within  the  anal 
area.  This  usual  positioning  of  the  swellings  is  shown  in  PI.  I,  A  and  B. 
A  few  swellings  occurred  within  the  anal  area,  and  these  were  frequently 
associated  with  a  minute  swelling  in  a  terminal  or  subterminal  position 
(PI.  I,  D).  In  certain  crosses  there  was  a  tendency  for  the  swellings  to 
be  more  distal  in  position  than  usual,  as  shown  in  PI.  I,  C.  Four  speci¬ 
mens  occurred  with  a  swelling  on  one  forewing  only. 

A  number  of  other  variants  of  the  forewing  occurred  with  some 
frequency  in  certain  crosses.  Some  of  these  are:  divergent  wings,  wings 
with  branching  anal  veins,  wings  shriveled  or  unfolded,  dwarfed  wings, 
and  wings  with  one  or  more  minute  blisters.  A  swelling,  filled  with 
hemolymph,  was  also  present  on  one  or  both  sides  of  the  pronotum  in 
some  balloon-wing  roaches,  but  this  character  does  not  appear  to  be 
confined  to  particular  crosses. 

On  the  hindwings  the  swellings  are  elongate,  usually  from  0.20- 
0.26  cm.  in  length  and  0.10-.12  cm.  in  width.  Sometimes  the  entire  wing, 
with  the  exception  of  the  anal  area,  is  inflated.  Hindwings  of  869  bal¬ 
loon-wing  roaches  were  examined.  Of  these,  818  had  swellings  on  one 
or  both  wings.  In  83%  of  the  specimens  the  swellings  occurred  within 
the  proximal  two-thirds  of  the  wing,  but  not  within  the  anal  area,  as 
shown  in  PI.  II,  A  and  B.  Swellings  within  the  anal  area  (PI.  II,  B)  or 
in  a  terminal  or  subterminal  position  (PI.  II,  C)  were  relatively  rare,  and 
confined  to  a  few  crosses.  There  were  103  roaches  with  one  normal¬ 
appearing  hindwing  (98 a* and  5  +  ),  and  approximately  half  this  number 
(51,  all  o»  )  which  had  no  swellings  on  either  hindwing. 

Several  variants  of  the  hindwing  occurred  which  may  be  of  genetic 
interest.  As  mentioned  in  the  preceeding  paragraph,  two  types  of  bubble¬ 
like  swellings  other  than  the  usual  balloon-wing  were  found.  In  addition, 
there  were  a  number  of  tiny,  blister-like  swellings,  a  broken  or  incom¬ 
plete  cubitus  (PI.  II,  C)  and  several  other  aberrations  of  the  venation. 

GENETICS 

The  balloon-wing  roaches  used  for  crossing  came  from  a  colony 
which  apparently  was  breeding  true  for  this  trait.  All  individuals  from 
the  colony  which  were  examined  showed  the  balloon  wing.  They  were 
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derived  from  those  few  mutants  originally  observed  in  a  normal-appear¬ 
ing  inbred  stock.  Those  roaches  involved  in  the  crosses  had  the  balloon  in 
the  usual  position  (PI.  I,  A  and  B)  with  a  few  in  PI.  I,  C.  The  wild 
type  (normal)  roaches  came  from  a  different  stock.  The  following  series 
of  crosses  were  made: 

Series  A.  —  Twenty-one  crosses  involving  balloon-wing  males  and  fe¬ 
males. 

Series  B.  —  Twelve  crosses  reciprocally  mating  balloon-wing  specimens 
with  normal  roaches. 

Series  C.  —  Twenty-four  crosses  reciprocally  mating  balloon-wing  speci¬ 
mens  and  normal-appearing  hybrids. 

Series  D.  —  Nineteen  crosses  involving  hybrid  males  and  females  ob¬ 
tained  from  Series  B. 

Series  E.  —  Nine  crosses  reciprocally  mating  normal-appearing  hybrids 
and  roaches  from  the  normal  stock. 

Result  of  Crosses 

The  results  obtained  from  the  five  series  of  crosses  listed  above  are 
presented  in  table  1.  The  reciprocal  matings  showed  no  significant  dif¬ 
ference  and  are  not  listed  separately  in  the  table.  The  results  indicate 
that  the  mutation  under  study  is  a  simple  Mendelian  recessive.  As  would 
be  expected,  balloon  x  balloon  (ser.  A)  produced  all  balloon  offspring; 
balloon  x  normal  (ser.  B)  produced  all  normal-appearing  offspring; 
balloon  x  hybrid  (ser.  C)  gave  results  which  closely  approximate  a 
ratio  of  one  balloon  to  one  normal;  hybrid  x  hybrid  (ser.  D)  produced 
offspring  in  a  ratio  of  three  normal  to  one  balloon;  and  hybrid  x  normal 
(ser.  E)  gave  all  normal  appearing  offspring.  The  results  from  series  C 
and  D  were  analyzed  statistically  by  the  chi  square  goodness  of  fit  test. 
It  was  found  that  the  chi  square  values  were  very  low.  This  indicates 
that  the  deviations  from  a  perfect  3:1  or  1:1  ratio  are  minor  and  are  well 
within  the  range  of  deviations  expected  by  chance  alone. 

The  four  specimens  mentioned  in  the  preceding  section  which  had 
swellings  on  only  one  forewing  are  believed  to  be  the  result  of  reduced 
expression  of  the  trait.  In  order  to  test  this  hypothesis  two  of  the  speci¬ 
mens  were  successfully  crossed,  one  with  a  normal  and  the  other  with  a 
balloon- wing  roach.  In  the  Fi  of  the  former,  all  offspring  were  of  normal 
appearance.  In  the  F2  there  were  190  normal-appearing  roaches  and 
65  roaches  with  swellings  on  both  forewings.  This  is  also  well  within 
the  expected  chance  deviations  from  a  3:1  ratio.  In  the  second  cross, 
there  were  26  males  and  32  females,  all  with  the  usual  balloon-wing 
swelling  on  both  forewings.  Thus  it  would  appear  that  in  these  two  indi- 
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viduals,  at  least,  the  absence  of  a  swelling  on  one  forewing  was  a 
phenotypic  variation. 

In  B.  germanica,  as  in  other  Orthoptera,  the  sex  determination 
mechanism  is  of  the  XO  type.  If  it  is  assumed  that  approximately  half 
of  the  sperm  carry  the  X  chromosome,  the  sex  ratio  should  be  1:1.  In 
the  present  study  an  opportunity  for  checking  this  ratio  was  present. 
The  results  gave  a  total  of  2132  males  and  2085  females  from  85  mat¬ 
ings.  This  is  well  within  the  expected  range  of  variation  for  a  1:1  ratio. 

DISCUSSION 

In  the  present  paper  the  character  “balloon-wing”  has  been  describ¬ 
ed  and  its  inheritance  elucidated.  The  associated  problem  of  how  the 
gene  responsible  for  this  trait  causes  the  swelling  to  develop  has  not  been 
investigated.  However,  the  writers  have  observed  the  development  of 
the  swelling  in  several  newly  emerged  adult  roaches  and  found  that  it 
first  shows  as  one,  two  or  three  tiny  blisters  which  appear  suddenly 
when  the  wing  is  partially  unfolded.  The  blister  or  blisters  enlarge 
quickly  to  form  a  swelling  of  the  usual  size  by  the  time  unfolding  is 
complete.  Normal-appearing  veins  are  present  before  the  appearance  of 
the  blisters  and  after  formation  of  the  swelling  are  still  discernable  on 
the  lower,  if  not  the  upper,  surface  of  the  swelling.  The  only  explana¬ 
tion  for  the  trait  which  can  be  offered  at  this  time  is  the  speculative  one 
that  perhaps  it  is  an  inheritable  weakness  of  the  middle  membrane  of 
the  wing. 

Along  these  same  lines,  it  is  interesting  to  note  that  the  character 
(balloon-wing)  does  not  always  express  itself  fully  in  roaches.  This  is 
probably  the  result  of  variable  expressivity  as  defined  by  Srb  and  Owen 
(1953:303).  Evidence  on  this  point  comes  from  the  two  crosses  reported 
in  the  previous  section  where  the  roaches  showed  a  swelling  on  only  one 
forewing.  In  both  cases  they  were  genetically  of  the  balloon  type.  If  it  is 
assumed  that  this  is  the  case,  then  a  calculation  of  the  per  cent  expression 
can  be  made.  For  this  purpose  it  will  be  considered  that  complete  ex¬ 
pression  prevails  in  all  those  specimens  in  which  a  balloon  appears  on 
both  wings  of  a  pair.  Accordingly,  full  expression  of  the  trait  in  fore¬ 
wings  is  greater  than  99  per  cent,  while  in  hindwings  it  is  approximately 
80  per  cent.  This  difference  in  expressivity  in  fore-  and  hind-wings  has 
not  been  studied,  but  it  may  be  genetically  controlled.  On  the  other 
hand,  there  is  an  explanation  which  does  not  require  any  genetic  involve¬ 
ment.  The  hindwings  normally  lie  in  a  folded  condition  between  the 
dorsal  body  wall  and  the  forewings.  It  is  possible  that  the  folds  and  the 
confined  position  confer  enough  structural  reinforcement  to  the  hind¬ 
wings  that  the  swellings  fail  to  appear  in  a  higher  percentage  of  hind¬ 
wings  than  forewings. 
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As  adults,  balloon-wing  Drosophila  lose  the  characteristic  swelling 
and  are  best  recognized  by  divergent  wings  and  an  extra  plexus  of  veins 
(Waddington,  1940:  128).  Balloon-wing  roaches  seldom  have  divergent 
wings  and  show  no  extra  plexus  of  veins.  The  swelling  on  the  wings 
are  the  only  constant  phenotypic  expression  of  the  mutation. 

SUMMARY 

A  new  mutation  has  been  described  for  the  German  cockroach, 
Blattella  germanica  (L.).  Based  on  analogy  with  similar  traits  in  other 
insects,  it  has  been  named  “balloon-wing.”  The  character  has  its  pheno¬ 
typic  expression  as  a  swelling  on  each  of  the  four  wings  of  this  insect. 
The  swellings  are  relatively  large,  and  on  the  forewings  usually  occur 
within  the  proximal  half  or  near  the  center  of  the  wing.  On  the  hind- 
wings  the  usual  expression  is  within  the  proximal  two-thirds  of  the  wing. 
In  a  small  percentage  of  cases  the  swellings  are  unusually  located.  A 
number  of  other  unusual  variants,  some  of  which  are  now  under  study, 
have  also  been  described. 

In  order  to  determine  the  mechanism  by  which  this  trait  is  inheri¬ 
ted,  a  series  of  crosses  involving  the  balloon-wing  stock  and  a  normal 
stock  have  been  made.  The  results  showed  that  the  swellings  are  the 
phenotypic  expression  of  a  simple  Mendelian  recessive  gene.  Four  speci¬ 
mens  (<1%  of  total)  were  found  with  a  swelling  on  only  one  fore¬ 
wing.  A  number  of  specimens  (approx.  20%),  usually  males,  were 
found  which  lacked  swellings  on  one  or  both  hindwings.  This  is  probably 
the  result  of  variable  expressivity,  and,  if  so,  the  per  cent  of  full  expres¬ 
sion  of  the  trait  is  considerably  lower  in  the  hindwings  than  in  the  fore¬ 
wings.  A  possible  explanation  for  this  difference  together  with  a  dis¬ 
cussion  of  several  other  points  is  included. 
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Table  1.  —  Results  from  a  series  of  crosses  involving  a  balloon-wing 
and  a  normal  stock  of  B.  germanica. 


Series 

No.  of 

normal-appearing 
offspring 
cf  °  total 

No. 

of  balloon-wing 
offspring 
a*  °  total 

A 

( balloon  x  balloon ) 

0 

0 

0 

456 

424 

880 

B 

(balloon  x  normal) 

417 

423 

840 

0 

0 

0 

C 

(balloon  x  hybrid) 

255 

247 

502 

279 

236 

515 

D 

(hybrid  x  hybrid) 

433 

440 

873 

140 

140 

280 

E 

(hybrid  x  normal) 

152 

175 

327 

0 

0 

0 
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Explanation  of  Plate  I 

Plate  I.  The  drawings  represent  left  forewings  of  B.  germanica,  and 
show  the  various  expressions  of  the  “balloon-wing”  mutation.  A  and  B, 
balloon  in  usual  location  on  wing;  C,  balloon  more  distal  than  usual; 
D  a  small  terminal  swelling  plus  a  balloon  in  the  anal  area. 
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Explanation  of  Plate  II 

Plate  II.  Similar  to  Plate  I.  Drawings  of  left  hindwings.  A,  balloon  in 
usual  position;  B.  balloon  in  the  usual  position  plus  a  balloon  in  the 
anal  area;  C,  balloon  more  distal  than  usual  plus  an  interrupted  cubitus 
vein. 
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THE  NEGATIVE  BINOMIAL  DISTRIBUTION: 

A  REVIEW  OF  PROPERTIES  AND  APPLICATIONS* 


John  J.  Bartko 
Virginia  Polytechnic  Institute 

1 .  Introduction 

The  negative  binomial  distribution  (NBD)  was  discussed  by  Pascal 
and  Fermat  (Todhunter,  1865).  However,  the  earliest  known  derivation 
and  publication  is  due  to  Montmort  in  1714.  A  particular  form  of  the 
negative  binomial,  the  inverse  binomial  sampling  form,  is  sometimes 
referred  to  as  the  Pascal  distribution. 

Many  forms  of  the  NBD  appear  throughout  the  literature.  The  most 
frequently  occurring  forms,  as  given  in  Table  I,  are  discussed  in  this 
paper. 


TABLE  I 


Forms  of  the  NBD 

Parameters 

Comment 

(D 

/X  +  r  -  1  >  r  x 
(  X  >P  1 

L  P 

Pascal  form;  p  +  q=l, 
r  is  an  interger. 

(2) 

/  n  -  lx  r  x 

(r  -  1>P  <1 

r>  P 

Pascal  form;  n=r+x. 

(3) 

/ x  +  k  -  lx  x  -x-k 

(  x  )P  q 

k,  p 

q=l+p;  k  not  neces¬ 
sarily  an  integer. 

The  NBD,  a  two  parameter  distribution,  has  the  same  sample 
space  (0,  1,  2,  .  .  .)  as  the  Poisson.  In  particular,  the  Pascal  form  (a 
term  reserved  specifically  for  the  case  where  r  is  an  integer)  is  a 
distribution  over  the  non-negative  integers. 

In  more  recent  years  NBD  applications  have  increased,  particularly 
in  health  and  accident  statistics,  binomial  sampling,  population  counts, 
data  of  contagious  distributions  and  communications. 

2.  Analytic  Properties 

Although  many  forms  of  the  NBD  are  utilized,  the  Pascal  form  is  of 
major  importance.  Its  primary  analytic  properties  follow. 
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The  characteristic  function  of  the  Pascal  form  is 

^x(t)  i  pr[(l  -  qelt)-r].  Accordingly,  the  mean  is 
rqp-1  and  the  variance  rqp-2 .  In  this  distribution  the  variance 
is  larger  than  the  mean.  From  the  form  of  $x(t)  it  is  easy  to  see 
that  the  sum  of  two  independent  negative  binomial  variates  x^  +  x2 
with  the  same  p  is  also  a  negative  binomial  variate. 

The  cumulant  generating  function  is 

^x(t)  =  /n  <J)x(t) 

=  -r  /n[  1  -  q(elt:  -  l)/p]  . 


From  the  definition  i.e.,  -  fi[x(x-l) . . .  (x-k+1)  ]  , 

the  factorial  moments  are  given  by  M-[kj  =  r(  r+1)  . . .  ( r+k-1)  p“kqk . 

A  very  interesting  and  convenient  property  of  the  Pascal  form  is  that 
its  c.d.f. 


(4) 


F(y) 


=  l  (: 

x=0 


+  r 
x 


r  x 
P  q 


=  Ip(r,  y+1) 


can  be  computed  bv  entering  tables  of  the  incomplete  Beta  function  as 
was  pointed  out  bv  Pearson  and  Fieller  (1933).  A  presentation  of  their 
result  follows. 

The  remainder  after  y  terms  of  the  Maclaurin  series  in  the  integral 
form  for  a  function  g(x)  is 

R  =  j  g^  (x-t)ty_1dt  . 

0 

Writing  the  Pascal  form  (1)  as  pr[(l-q)“r]  and  taking 
g(q)  =  (l-q)“r  we  have 

R  =  T$frfy)  f  (i-q^t)-r'yty-1dt  . 

J0 

Let  (l-u)  =  (l-q)/(l-q+t) ;  (5)  with  the  factor  pr  included 

then  becomes 

0 
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and  finally  R  =  I^(y,  r)  , 

where  Ix(a,  b)  is  the 

standard  notation  for  the  incomplete  Beta  function 
x 

ta-1(l-t)b_1dt/&(a,b). 

0 

Therefore, 

^  x  +  r  -  IjpT  x.  _  j  (y^r) 
x=y  x  q 

and 

(6)  Fly)  =  1  -  Iq(y+1,  r)  =  Ip(r,  y+1)  . 

It  is  interesting  to  note  that  the  use  of  incomplete  Beta  tables  for 
the  evaluation  of  this  c.d.f.  has  been  rediscovered  by  Patil  (1960). 

3.  Limiting  Forms 

Three  limiting  forms  of  the  NBD  are  discussed  in  this  section. 

3.1  The  Poisson  Distribution. 

Assume  the  Pascal  form  ( 1 )  whose  moment  generating  function 
(m.g.f.)  is  pr/(l-qe9)r.  Let  X.  be  fixed  such  that  q  =  X./r, 
and  let  p  -*  1,  q  -»  0  then 

(7)  (-P  Q)r  =  , 

r  -#•  oo  1-qe 

the  m.g.f.  of  the  Poisson  distribution,  e_:Vx/x!,  where  x  is  the 
number  of  failures  before  the  desired  number  of  successes  r. 

3.2  The  Gamma  Distribution. 

Consider  the  Pascal  form  (2)  giving  the  distribution  of  sample  size 

n  as 

f ( n)  =  - ^n~1^ -  prqn-r  ,  n  =  r,  r+1,  ... 

(n-r) ! (r-1) ! 

Suppose  each  trial  takes  time  At  and  that  the  probability  of  success  p 
is  proportional  to  At ,  i .  e . ,  p  =  AAt .  then  time  t  =  nAt  required 
to  achieve  r  successes  has  the  d.f.  f(t)  given  by 
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f(t)at  =  -5 r — 


(tt-  -  1)  !• 


(Mt)r(l  -  Mt)At 


or 


f(t)  =  (t  -At)...(t  -  r-lAt )^r( 1-X.At ) 
(r  -  1) ! 


At 


Let  At  •+  0;  then  f(t)  becomes  the  Gamma  density 


(8) 


•x  r  r-1  -\t 

f(t)  =  - — - — - — 

r(r) 


which  may  be  interpreted  as  the  “holding  time”  before  the  rth  occurrence 
of  a  rare  event. 

3.3  The  Logarithmic  Distribution. 

In  many  applications  of  the  NBD 


(9) 


(k+i-Dlp1 

(k-1) !i ! (l+p)k+i 


i  0,1,2,... 


the  number  i  observed  in  any  sample  may  have  all  integral  values 
including  zero.  However,  in  some  cases  the  value  i  =  0  will  not  be 
observed.  For  example  in  the  application  of  the  NBD  to  the  number  of 
representatives  oi  different  species  of  butterflies  obtained  in  a  collection, 
only  frequencies  of  numbers  greater  than  zero  will  be  observable,  since 
by  itself  the  collection  gives  no  indication  of  the  number  of  species 
which  are  not  found  in  it. 


The  following  development  of  the  logarithmic  distribution  is  due  to 
Fisher  (1943).  In  nature  the  abundance  of  different  species  of  the  same 
group  generally  varies  very  greatly.  Fisher  found  in  a  study  of  Malavan 
butterflies  that  the  NBD  which  often  fits  such  data  well  has  a  value  of 
k  so  small  as  to  be  almost  indistinguishable  from  zero.  That  it  is  not 
really  zero  for  collections  of  wild  species  follows  from  the  fact  that  the 
total  number  of  species  and  therefore  the  total  number  not  included 
in  the  collection  is  really  finite.  The  real  situation  in  which  a  large 
number  of  species  are  so  rare  that  their  chance  of  inclusion  is  small,  is  well 
represented  by  the  limiting  form  of  the  NBD  when  k  tends  to  zero. 

The  limiting  value  k  =  0  cannot  occur  in  cases  where  the  frequency 
at  zero  is  observable,  for  the  distribution  would  then  consist  wholly  of 
such  cases.  If,  however,  we  put  k  =  0  in  (9),  write  x  for  p/(l+p),  so 
that  x  is  a  positive  number  less  than  unity  varying  with  the  size  of  the 
sample  and  replace  the  constant  factor  (k-1)!  in  the  denominator  by  a 
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new  constant  factor  a  in  the  numerator,  we  have  an  expression  for  the 
expected  number  of  species  with  i  individuals,  where  i  now  cannot  be 
zero, 

ax1/ i . 

The  total  number  of  species  expected  is  consequently 

00  -i  „ 

(10)  E  ax  /i  =  -a  /n(l-x) 
i=l 

the  logarithmic  distribution.  The  total  number  of  individuals  expected  is 
00  • 

(11)  E  ax1  =  ax/(l-a). 
i=l 

Thus  the  two  relationships  (10)  and  (11)  enable  the  series  to  be  fitted 
to  observational  data,  for  if  S  is  the  number  of  species  observed,  and  N 
the  number  of  individuals  the  two  equations 

(12)  S  =  -  a  jfn(l  -  x)  ,  N  =  ax/(l  -  x)  , 

are  sufficient  to  determine  the  values  of  a  and  x.  Tables  which  facilitate 
the  solution  of  these  equations  and  an  example  can  also  be  found  in 
Fisher  (1943). 


4.  Truncated  NBD 

In  fitting  data  where  the  zero  class  is  not  observed  it  may  be  the 
case  that  a  truncated  NBD  gives  a  better  fit  than  the  logarithmic 
distribution.  Sampford  (1955)  considers  the  truncated  NBD  and  gives 
procedures  and  tables  to  obtain  iteratively  either  the  moment  or  maximum 
likelihood  estimates  of  the  parameters.  Included  also  is  a  discussion  of  the 
efficiency  of  these  procedures. 

Simplified  methods  for  the  estimation  of  the  parameters  of  the 
truncated  distribution  by  a  modified  method  of  moments  and  by  a 
modification  of  maximum  likelihood  procedures,  due  to  Brass  (1958) 
are  discussed  below.  Brass  also  discusses  the  efficiency  of  these  two 
procedures. 

Consider  the  NBD  of  the  form  (3) 


(13) 


Pi 


(k+i-1) Ip1 


(k-1) ! i ! ( 1+p ) 
so  that  P(o)  =  l/(l+p)k. 
bilities  (13)  must  be  divided  by 


k+i 


,  i  =  0,  1,  2,  . . . 

To  obtain  the  corresponding  proba- 
[1  -  P(0)];  writing  u>=i/(i+p) 
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and  T)  =  1  —  oo  it  follows  that 

pTu)  =^Hk±iziiinL  , 

T  (l-ook)(k-l)  !i! 


i  =  1,  2,  ... 


where  the  subscript  T  denotes  the  truncated  distribution. 

By  considering  the  proportion  in  the  first  class  of  the  truncated 
distribution  and  its  mean  and  variance  found  by  use  of  factorial  moments 
Brass  develops  simple  estimates  for  oo,  k  and  also  for  M,  the  mean  of  the 
complete  NBD.  The  value  M  enters  in  his  discussion  of  efficiency.  Aside 
from  M,  the  estimates  are 


=  -m(l  -  n-L/n)/s: 


(14) 


k  = 


OQm  -  n^/n 

~  } 

1-00 


where  the  bars  denote  the  estimates.  m  is  the  number  of  sample 
observations  with  measurement  i  and  n  is  the  total  sample  size,  m  and 
s2  are  the  unbiased  estimates  of  |i  and  a2,  that  is 


oo  in. 
m  =  £  — - 

i=l  n 


oo  n .  (  i-m)  2 
E  1 


i=l 


n  -  1 


Considering  the  log  likelihood 


Jfn  L  -  Z  n.jfnC^^DKl-")1] 
1  (l-cok)  (k-1)  !i! 


i=l 

Brass  obtains  the  following  modified  maximum  likelihood  expressions, 
k  +  n^/n 


(15) 


k  +  m 


(16) 


m(k  +  n-,/n)  k  +  n-,/n 

= - V-  /n  (— - ) 

k(m  -  n-j^/n)  k  +  m 


1 

n 


R  _ 

Z(k  +  j 


1) 


i=j 


n. 

l 


0 


where  R  is  the  largest  observed  value  of  i.  The  estimate  from  (14)  gives 
a  convenient  first  entry  for  k  in  the  iterative  procedure  (16).  co  then 
follows  immediately  from  (15). 

The  above  methods  of  fitting  are  illustrated  on  the  data  below  also 
from  Brass.  The  observations  are  the  number  of  children  ever  born  to  a 
sample  of  mothers  over  forty  years  of  age. 
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No.  of  Children  Per  Mother 

i 

1  23456789  10  11  12 

Total 

No.  of 

Mothers,  n. 

1 

49  56  73  41  43  23  18  18  7  7  3  2 

340 

m  =  3.9912,  s2  =  5.9734,  n±/n  =  49/340  =  0.1441- 


Then  from  (14)  co  ;=  0.572  and  k  =  5. 00. 

For  the  maximum  likelihood  solutions  write  (16)  as 


m-n-jyn 


(1  + 


nn/n 


)/n(l  - 


/n 


k+m 


)  +  i  I  (k+j-1) 
j=l 


-1 


h  n. 

i-d  1 


Then 

i  -  4>* 


_  mn,  / n 

m  V _ 

k ( k+m )  "  k2(m-n^/ n. 


/n(i_ 


m-n 

k+m 


n 

-) 


-in  9  n 

-  2  (k+j-1)  £  n  . 

J-l  i=j 


The  iterative  procedure  is  to  calculate  (|)  and  (J)T  for  some  value  of  k; 
a  second  value  of  k  which  will  make  ())  closer  to  zero  is  then  found  by 
the  aid  of  (j)'  on  the  assumption  that  the  relation  between  (J)  and  k  is 
roughly  linear  over  a  short  range.  Further  improved  approximations  are 
obtained  by  linear  interpolation  from  the  values  of  $  for  the  preceding 
two  estimates  of  k. 


k 

4> 

4>f 

5.00 

-0.00064 

-0.00258 

4.75 

+  0.00008 

- 

4.78 

0.00000 

- 

The  first  approximation  to  k  is  taken  as  5.00  from  the  estimate  by  (14) 
The  second  approximation  to  k  is 


5.00 


-0.00064 

-0.00258 


=  4.75. 


One  linear  interpolation  between  5.00  and  4.75  gives  4.78.  From 
k  =  4.78,  can  be  found  directly  from  (15)  giving  Za  =  0.561.  The 
estimates  again  from  (14)  are  k  =  5.00  and  00  =  0.572. 

The  estimates  by  the  two  procedures  differ  little.  Brass  gives  the 
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efficiencies  for  both  methods  at  about  95%.  These  efficiency  values  are 
read  from  tables  which  appear  in  his  paper.  The  tables  are  based  on 
estimates  of  k  and  M.  In  general  the  modified  moments  method  of 
estimation  is  more  valuable  in  that  slight  labor  is  involved  and  efficiency  is 
high  in  many  cases  met  in  practice. 

5.  Occurrence 

Inverse  binomial  sampling  is  of  primary  importance  in  the  considera¬ 
tion  of  the  NBD.  It  occurs  in  the  following  general  type  of  situation.  If 
a  proportion  p  of  individuals  in  a  population  possesses  a  certain 
characteristic,  the  number  of  observations  in  excess  of  r  which  must  be 
taken  to  obtain  exactly  r  individuals  with  the  characteristic  has  the  Pascal 
form  of  the  NBD. 

Greenwood  and  Yule  (1920)  study  the  NBD  in  connection  with  a 
Poisson  distribution  with  varying  mean. 

Yule  (1924),  Furry  (1939)  and  Kendall  (1949)  consider  the  dis¬ 
tribution  in  connection  with  birth-death  processes.  Yule  is  concerned  with 
the  mathematical  theory  of  evolution;  Furry  with  cosmic  radiation  and 
Kendall  with  the  process  in  general. 

McKendrick  (1914)  and  Quenouille.  (1949)  discuss  contagious  dis¬ 
tributions.  Irwin  (1941)  and  Luders  (1934)  are  concerned  with  the 
NBD  as  it  occurs  in  multiple  accident  theory. 

The  first  of  these  situations,  inverse  binomial  sampling,  is  the 
simplest  mathematically  and  is  the  only  one  where  the  mathematical  model 
is  likely  to  hold  exactly  in  practice. 

6.  Computational  Procedures  on  the  Pascal  Form 
6.1  Point  Estimate  of  p. 

Occasions  may  arise  in  which  the  frequency  of  an  attribute  p  is 
unknown.  For  the  Pascal  form  the  estimation  of  p  from  the  sample  data 
is  a  simple  matter. 

Haldane  (1943)  considers  the  NBD  (2)  and  derives  an  unbiased 
estimate  of  p,  viz. 

(17)  r-1 

n-1 

where  n  is  the  total  number  of  items  sampled  in  order  to  recover  the 
predetermined  number  r.  Haldane  also  develops  an  expression  for  the 
variance  of  this  estimate.  It  is,  however,  approximate.  Finney  (1949) 
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derives  an  exact  expression  for  the  variance  and  also  obtains  an  unbiased 
estimate  of  the  variance. 

6.2  Interval  Estimate  of  p. 

It  may  not  have  been  generally  realized  that  methods  and  tables  for 
determining  exact  confidence  limits  on  p  in  direct  or  binomial  sampling 
may  be  adapted  very  easily  to  inverse  binomial  sampling.  The  proof 
which  follows  is  due  to  Finney  (1949). 

If  pu  is  the  upper  confidence  limit  for  a  value  of  p  estimated  from 
inverse  sampling,  by  definition  it  is  the  value  of  p  for  which  the  proba¬ 
bility  of  obtaining  the  quota  r  in  n  or  more  trials  is  exactly  Ja.  This 
value  pu  is  the  same  value  of  p  which  would  give  r-1  or  less  with  the 
attribute  in  a  sample  of  predetermined  size  n-1  in  direct  sampling.  Like¬ 
wise  pL,  the  corresponding  lower  limit,  is  the  value  of  p  which  would 
give  the  quota  in  n  trials  or  less  with  a  probability  of  ,  and  is,  there¬ 
fore,  the  same  as  the  value  of  p  which  in  direct  sampling  would  give  r  or 
more  with  the  attribute  in  a  sample  of  size  n.  Thus,  the  rules  for  adapt¬ 
ing  Table  41  of  Bometrika  Tables  for  Statisticians  which  gives  confidence 
limits  for  the  direct  binomial  can  be  stated: 

(i)  For  the  upper  confidence  limit  on  p  enter  the  table  with 
c  ==  r-1,  n  =  n-1. 

(ii)  For  the  lower  confidence  limit  enter  the  table  with  c  =  r, 
n  =  n. 

As  an  example  of  the  use  of  these  rules  consider  the  following 
situation.  An  experimenter  has  concluded  an  inverse  sampling  plan  in 
which  he  obtained  four  desired  attributes  in  ten  trials.  Here  r  =  4, 
n  =  10.  To  find  the  upper  confidence  limit  on  p  enter  the  table  with 
c  =  3,  n  =  9.  For  the  lower  confidence  limit  enter  with  c  =  4,  n  =  10. 
Consequently,  the  two  sided  95%  confidence  interval  for  p  is 
0.11  <  p  <  0.70.  The  point  estimate  of  p  from  (17)  is  1/3. 

6.3  An  Example  of  Inverse  Binomial  Sampling. 

Chapman  (1952)  considers  an  inverse  binomial  sampling  plan  in 
which  it  is  desired  to  estimate  the  size  N  of  a  population.  In  Chapman’s 
plan  a  number  t  of  individuals  are  taken  from  this  population.  They  are 
tagged  and  returned.  In  the  second  sample,  n  individuals  are  sampled  in 
order  to  recover  a  predetermined  number  r  of  the  tagged  ones. 

If  sampling  is  with  replacement  P(n)  the  probability  of  having  to 
sample  n  individuals  to  obtain  r  tagged  ones  is  given  by 


p<n)  =  -  H)n_r 
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where  r  is  the  fixed  parameter  and  n  the  random  variable.  £(n)  =  rN/t 
and  nt/r  is  an  unbiased  estimate  of  N.  Chapman  considers  a  normal 

approximation  and  a  x2  approximation  for  confidence  limits  on  N. 
However,  a  more  expedient  method  is  to  use  the  rules  stated  in  (6.2)  and 
Biometrika  Table  41.  Let  t/N  =  p,  replacing  N  by  nt/r,  p  becomes  r/n. 

Then  pn  and  pL  can  be  determined  such  that  p(Pl  -  P  5  p^j)  >  1-cl. 

On  replacing  p  by  t/N  the  confidence  statement 
P^/p-pf  -  N  -  t^pL^  -  1-a  obtained. 


7.  Estimation  of  Parameters 


The  estimation  of  parameters  of  other  forms  of  the  NBD  is  not  as 
straightforward  as  it  is  in  the  case  of  the  Pascal  form.  Parameter  estima¬ 
tion  in  these  cases  is  usually  accomplished  by  the  method  of  moments. 
However,  when  this  method  is  not  satisfactorily  efficient,  the  method 
of  maximum  likelihood  is  commonly  used. 


Consider  first  the  method  of  moments.  Fisher  (1941)  discusses  the 
NBD  (3) 

f(x)  = 


(k+x-1) !p" 


(k-1) !x! ( 1+p) 


k+x 


and  finds  the  estimates  of  p  and  k  by  the  method  of  the  moments.  These 
estimates  are 

P*  =  (s2  -  x)/x 


(18)  and 

k*  =  x2/(s2  -  x)  . 

Fisher  also  compares  the  determinants  of  the  covariance  matrices  for 
large-sample  estimates  found  by  moments  and  by  maximum  likelihood 
obtaining  an  expression  for  the  reciprocal  of  the  efficiency.  This  further 
yields  an  expression,  viz. 

(19)  (1  +  1/p)  (k  +  2) 

which  indicates  when  moment  estimates  are  satisfactorily  efficient,  i.e., 
>  90%.  This  is  the  case  when  (19)  exceeds  20  using  the  values  of  p 
and  k  found  from  (18).  If  the  guide  (19)  indicates  that  moment 
estimates  are  not  satisfactory  then  another  method  of  estimation  must  be 
employed. 

Estimation  by  maximum  likelihood  is  efficient.  Haldane  (1941) 
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considers  the  estimates  of  p  and  k  by  maximum  likelihood.  The  estimate 
of  p  is 

(20)  P  =  x/k 


where  k  is  obtained  iteratively  from  the  following  expression 


(21)  N[jTn(k+x)-jfn(k)]  = 


n1+n2+---+nR 


n2+...+nR 

k+1 


nr 


+  .  .  .+ 


k+R-1 


ru  is  the  observed  frequency  of  xi  ,  R  is  the  the  maximum  value  of  xi 


such  that  i 

of  observations  and 


0,  1,  R. 

-  1  R 
X  -  N 


R 

N  =  2  n.  , 

i=0  1 


the  total  number 


As  an  illustration  of  the  two  methods  of  estimation,  consider  the 
following  data  from  Fisher  (1941).  The  table  gives  a  sample  of  sheep 
classified  according  to  the  number  of  ticks  found  on  each. 


i 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

number  of 
ticks  x- 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

number  of 
sheep  rq 

7 

9 

8 

13 

8 

5 

4 

3 

0 

1 

2 

2 

II 

02 

O 

x  =  3.25  s2  =  5.9195. 

From  (18)  k  =  3.75  and  p  ==  0.821.  With  these  values  (19)  is 
13.2  so  that  maximum  likelihood  estimates  are  the  desirable  ones.  By 
(21)  k  =  3.75.  p  =  0.867  follows  directly  from  (20). 

8.  Hypothesis  Testing 

Barnard’s  (1946)  method  for  testing  the  hypothesis  Hq:  p1  =  p2  in 
inverse  binomial  sampling  is  presented  below. 

Consider  the  situation  in  which  an  instrument  has  been  modified  in  the 
hope  of  removing  a  cause  of  failure.  The  sampling  procedure  consists  in 
arranging  to  try  the  modified  one,  under  similar  conditions,  and  to  con- 
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tinue  to  try  each  until  a  predetermined  number  rx  of  failures  of  the 
modified  instrument  and  a  predetermined  number  r2  of  failues  of  the 
unmodified  instrument  are  observed  (inverse  binomial  sampling).  Then 
if  p!  refers  to  the  modified  instrument  and  p2  to  the  unmodified  one,  and 
nx,  n2  are  the  number  of  observations  made  in  the  two  cases,  the  proba¬ 
bility  of  the  pair  (nx,  n2)  is 


n-,  -1 


1  J 
ri)pi 


n, -r,  n~-l  r0 
.  1 


n2~r2 


If  pi  =  p2  =  p,  the  probability  of  obtaining  some  pair  of  results  having 
the  same  total  nx  +  n2  =  N  is 


1  .  rl+r2  N-(rl+r2> 

)P  9 


Then  the  relative  probability  on  the  null  hypothesis  px  =  p2  of  obtaining 
the  pair  (nx,  n2)  out  of  all  results  with  he  same  total  N  is 


(22) 


LtMr-j+i^)  ]!  (r-j+r^l)  !(n1-l)  1(^-1)  ! 
( N-l) I (r1— 1) ! ( r2“l) 1 ^ni"ri^ ! ^n2"r2^ ! 


which  is  independent  of  p.  Hence,  if  we  sum  over  more  extreme  pairs 
having  the  same  N,  rx  and  r2  we  shall  obtain  the  significance  level  of  our 
result. 


For  example,  with  rx  =  1,  r2>  =  2  if  we  find  nx  =  10,  n2=4  we  have 
the  significance  level 


2 ! 11 !  ,  9 ! 3 !  ,  10 !2 ! 
13 roil!  ' 9T2T  1011! 


+ 


11!1I 

1110! 


)  =  1/13 


which  is  non-significant  at  the  5%  level. 

Barnard  also  develops  an  approximation  procedure  when  px  and  p2 
are  small  and  n!  and  n2  are  large.  The  rule  is  to  calculate 
R  =  [n-[-T^/2)T2/{r\2-r2/2)r1Sind  to  enter  the  F  tables  with  this  value  of 
R  for  (2r1,  2r2)  degrees  of  freedom. 

For  the  case  of  a  2xjf  trial,  provided  px,  p2  .  .  .,  Vf  are  all 

small,  Bartlett’s  or  Hartley’s  test  can  be  used.  Barnard  did  not  elaborate 
this  point.  It  seems  worth-while,  therefore,  to  examine  how  this  situation 
can  be  referred  to  Bartlett’s  test. 


1/2  )A- 


XIr. 


Note  that  2(n 


and 
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Identify  -  a2 


and 


(nl  -  T-J2)/Txn  . 

(rq  -  rq/2)/iq  =  s|r  i.e.,  s2  with  2r1  degrees  of  freedom. 

Bartlett’s  test,  M  =  (N*  -  jf)jfn  s2  -  S[(n!-l)/n  s2]  where  n»  -  total 

number  of  observation,  n!  =  number  of  observations  in  the  ith  sample, 
Jf  =  number  of  samples,  then  becomes  in  our  notation 


(23)  M  =  2Erq/n  s2  -  2Zrijfn  s? 

where  s?  =  (n±  -  r±/2)/r±,  i  =  1,  2,  . . . ,  jf 

and  s*  =  Eire  -  r±/2)/Zr±  . 


To  perform  the  test  refer  M  to 

If  ri  =  r2  =  *  *  *  =  r/  perform  Hartley’s  test,  max  s2/min  s2,  and 
refer  to  Biometrika  Table  31. 

9.  x2  -Tests  for  Contingency  Tables  of  Negative 
Binomial  Types 

Steyn  (1959)  discusses  the  problem  of  the  analysis  of  a  contingency 
table  where  the  frequency  of  one  of  the  cells  is  predetermined.  A 
summary  of  his  work  on  this  problem  follows. 

Suppose  sampling  with  replacement  is  performed  in  a  population 
consisting  of  categories  FE,  FE,  FE,  and  FE  in  proportions 
Pll'  Pol  ’  P]_o  and  Poo  resPectivety*  Individuals  are  drawn  at  random 
until  r  FE’s  are  found.  Obviously  the  sample  size  is  a  random  variable. 

The  following  table  shows  the  observed  frequencies  for  the  negative 
fourfold  table.  This  table  is  obtained  when  two  events  E  and  F  can  occur 
or  fail  either  alone  or  together  and  sampling  is  continued  until  the  rth 
joint  failure  has  occurred,  thereby  fixing  a  priori  the  frequency  in  one 
cell  but  not  the  total  in  the  sample. 


1961] 


Bartko— The  Negative  Binomial 
Frequencies 


31 
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E 

F 

r 

f01 

r  +  f01 

F 

0 

1 — 1 

fu 

fio  +  fii  -  y 

r  + 

f  10 

f01  +  fll  =  x 

Steyn  develops  the  expression  for  the  test  on  this  table.  The 
expression  with  one  degree  of  freedom  is 


(24) 


X2 


_ 2 _ 


01 


-*•#> 


r(&) 


(?) 


1 

'  *  ) 

kqqf 


where  (p?/qT)  =  [-(r-x+y)  +  /(r+x+y)‘i  -  4xy)]/2r  , 

(p/q)  =  [-(r+x-y)  +  /(r+x+y)2  -  4xy  ]/2r 

and  (i/qq»)  =  i  +  (pT/qf)  +  (p/q)  +  (p/q) (pT/qf ) • 


The  extension  to  the  case  where  there  are  h  possibilities  for  the  event 
E  and  jf  possibilities  for  F  is  also  discussed. 

Summary 

This  paper  is  a  review  of  the  major  literature  dealing  with  the 
negative  binomial  distribution  (NBD).  An  account  of  the  forms  of  the 
distribution  is  presented.  This  includes  limiting  and  truncated  forms. 
Parameter  estimation  is  treated  in  a  few  specific  cases.  The  problem  of 
testing  p-|_  =  P2>  also  that  of  testing  Pi  =  P2  =  *“  =  pjf  along 

with  a  discussion  on  contingency  tables  of  negative  binomial  types  is 
presented  in  some  detail. 

The  inverse  binomial  sampling  form  is  emphasized  throughout  the 
paper.  In  the  discussion  of  this  form  several  extremelv  useful  relations  are 
discussed.  Among  these  is  a  discussion  of  the  use  of  the  Tables  of  the 
Incomplete  Beta  Function  for  the  evaluation  of  the  cumulative  distribu¬ 
tion  function  (c.d.f)  Also  discussed  is  the  method  of  adapting  Table  41 
of  Biometrika  Tables  for  Satisticians,  giving  confidence  limits  on  p  in  the 
positive  binomial,  to  the  inverse  binomial  sampling  case. 
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McCOMAS 

RESEARCH 

CENTER 


to  provide  for  growing 
basic  research  program 


The  opening  of  our  new  research 
center  this  month  will  provide  us 
with  an  opportunity  for  even 
more  extensive  basic  research 
into  the  physical  sciences. 

The  new  laboratories  and  in¬ 
strumentation  will  facilitate 
comprehensive  studies  into  the 
basic  nature  of  tobacco ,  smoke 
and  smoke  flavor  as  well  as  in¬ 
tensify  our  corporate  research 


into  packaging  and  adhesive  ma¬ 
terials. 

We  are  looking  forward  to  the 
achievements  to  come  from  this 
new  laboratory  complex.  We 
sincerely  believe  our  greater  em¬ 
phasis  on  fundamental  research 
will  be  an  important  factor  in 
the  continued  growth  of  both 
our  company  and  the  industries 
we  serve. 
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NEWS  AND  NOTES 


( Editors  Notes  —  News  contributions  should  be  sent  to  the  person  whose 
name  appears  at  the  end  of  the  appropriate  sections. ) 

ENGINEERING  SECTION 

Stricken  by  heart  attacks  during  the  Fall,  both  Professor  Fred 
Bull  of  the  chemical  engineering  department  at  Virginia  Polytechnic 
Institute  and  Mr.  Tilton  E.  Shelburne  of  the  Virginia  Council  for 
Highway  Investigation  and  Research  are  reported  well  on  the  road  to 
recovery  and  expect  to  resume  their  duties  early  in  1961. 

Dean  Lawrence  R.  Quarles  of  the  School  of  Engineering,  University 
of  Virginia,  attended  a  conference  on  University  Reactors  at  Gatlinburg, 
Tennessee,  sponsored  by  the  Oak  Ridge  Institute  of  Nuclear  Studies 
and  presented  a  paper  “Construction  of  the  University  Reactor.”  In 
October  Dean  Quarles  attended  a  meeting  of  the  ANS-ASEE  Committee 
on  Objective  Criteria  for  Nuclear  Engineering  Education  in  Chicago  and 
also  a  meeting  of  the  Region  IV  Executive  Committee  of  the  AIEE  in 
Charlotte,  North  Carolina.,  Professor  Orville  R.  Harris  of  the  electrical 
engineering  department  of  the  University  of  Virginia  attended  a  sympo¬ 
sium  in  New  York  in  November  on  Magnetism  and  Magnetic  Materials. 
In  December  Dr.  Harris  attended  the  annual  meeting  of  the  American 
Rocket  Society  in  Washington. 

Professor  R.  E.  L.  Gildea  spent  a  portion  of  the  summer  of  1960  on 
active  duty  with  the  U.  S.  Public  Health  Service  during  which  time  he 
helped  organize  a  regional  meeting  on  Radiological  Health  to  be  spon¬ 
sored  during  the  spring  by  the  University  of  Virginia,  the  Virginia  De¬ 
partment  of  Health,  and  the  U.  S.  Public  Health  Service.  Major  D.  M. 
Crim  of  the  civil  engineering  department  at  Virginia  Military  Institute 
is  on  leave  of  absence  during  the  first  semester  of  the  1960-’61  academic 
year  and  is  studying  for  a  graduate  degree  in  civil  engineering  at  the 
University  of  Virginia. 

Professors  Otis  L.  Updike  and  Robert  M.  Hubbard  of  the  chemical 
engineering  department,  University  of  Virginia,  attended  the  annual 
meeting  of  the  A.  I.  Ch.  E.  in  Washington  in  December.  Professor 
Updike  was  co-chairman  for  two  sessions  of  the  symposium  on  “The 
Unsteady  State  as  an  Analysis  Tool  in  Chemical  Engineering.”  Professor 
Hubbard  is  chairman  of  two  subcommittees  of  the  A.I.Ch.E.  Chemical 
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PHIPPS  &  BIRD 

PIPETTE  DRYER 

In  the  laboratory  this  large 

capacity  Pipette  Dryer  is  special¬ 
ly  recommended  for  thorough  dry¬ 
ing  pipettes  of  any  size  up  to 
375  mm.  in  length. 

Cylinder  is  mounted  on  vented 
base  with  specially  designed  heat¬ 
er  placed  in  lower  end  of  cylin¬ 
der.  The  instrument  is  designed 
to  take  the  standard  rack  of 
6-in.  pipette  washer.  Rack  is  tak¬ 
en  out  of  washer  drained  for 
about  a  minute  and  placed  right 
into  the  dryer.  Drying  time  varies 
according  to  size  and  load  of 
pipettes.  But  this  method  of  dry¬ 
ing  means  less  handling  by  the 
technician  thus  the  hazard  of 
breakage  is  greatly  reduced. 


Cat.  No.  71-281-1  Cat.  No.  71-726 
$45.00  $10.00 


Without  obligation  send  us  com¬ 
plete  information  on  your  Pipette 
Dryer. 

Name  . . 

.  Title  . 

Company  . 

Address  . 

City  . 

.  State  . . 

The  Pipette  Dryer  comes  com¬ 
plete  with  cord  and  plug  but 
without  rack.  The  instrument 
operates  on  115  volts  AC  or  DC. 
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Engineering  Education  Projects  Committee.  A  paper  by  Professor  Hubbard 
‘‘Teaching  with  a  Pneumatic  Analog”  was  published  in  the  August  and 
September  issues  of  the  ISA  Journal. 

Professor  Bull  has  announced  that  Mr.  R.  H.  Pusey  has  assumed  a 
teaching  position  on  the  staff  of  his  department  at  VPI  taking  the  place 
of  Mr.  G.  F.  Meehan  who  recently  resigned. 

Professor  Nelson  F.  Murphy  of  the  chemical  engineering  depart¬ 
ment  at  VPI  is  a  coauthor  with  Professor  F.  W.  Bull  of  a  bulletin  from 
the  Virginia  Engineering  Experiment  Station  entitled  “Decolorization  of 
Waste  from  Caustic  Wash  of  Chlorine-Bleached  Kraft  Pulp.”  Professor 
Murphy  attended  a  meeting  of  the  Committee  on  Education  of  the 
American  Electroplaters’  Society  in  Hartford,  Connecticut  in  February. 
In  November  he  spoke  to  the  Blue  Ridge  Branch  of  the  American 
Electroplaters’  Society  in  Roanoke  on  the  subject  “How  Acid  Inhibitors 
Work/’ 

Professor  Robert  Truitt  of  Virginia  Polytechnic  Institute  has  an¬ 
nounced  that  the  name  of  his  Department  of  Aeronautical  Engineering 
has  been  changed  to  the  Department  of  Aerospace  Engineering. 

Mr.  Tilton  E.  Shelburne  of  the  Virginia  Council  for  Highway  Investi¬ 
gation  and  Research  attended  a  meeting  in  Chicago  in  February  of 
the  Data  Analysis  Panel  on  the  AASHO  Road  Test. 

Dean  Lawrence  R.  Quarles  of  the  Engineering  School  of  the  Uni¬ 
versity  of  Virginia  attended  a  meeting  of  the  American  Society  for 
Engineering  Education- American  Nuclear  Society  Committee  on  Ob¬ 
jective  Criteria  for  Nuclear  Education  in  Dallas,  Texas  in  January.  Also 
during  that  month  Dean  Quarles  attended  a  national  meeting  of  the 
American  Institute  of  Electrical  Engineers  in  New  York.  In  February  he 
attended  a  Conference  on  Higher  Adult  Education  sponsored  by  the 
Southern  Regional  Education  Board.  He  also  met  with  the  Regional 
Advisory  Council  on  Nuclear  Energy  in  Atlanta,  Georgia  in  February. 

Professors  Otis  L.  Updike  and  Robert  M.  Hubbard  of  the  chemical 
engineering  department  of  the  University  of  Virginia  attended  the 
annual  meeting  of  the  American  Institute  of  Chemical  Engineers  in 
Washington  in  December.  Dr.  Updike  was  co-chairman  of  two  sessions 
of  a  symposium  on  the  subject  “The  Unsteady  State  as  a  Tool  for  Analy¬ 
sis  in  Chemical  Engineering.”  Professor  Updike  has  recently  received  a 
subcontract  from  the  Oak  Ridge  National  Laboratory  for  research  on 
the  dynamics  of  solvent  extraction  cascades. 
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Drs.  D.  G  Cochran  and  J.  M.  Grayson,  V.  P.  I.  Entomologists,  have 
received  a  grant  of  $3,000  for  studies  in  cross  resistance  and  control  of 
German  cockroaches  resistant  to  organic  phosphate  insecticides.  The 
research  will  be  in  cooperation  with  the  National  Pest  Control  Associa¬ 
tion.  Both  investigators  were  speakers  at  recent  Regional  Pest  Control 
Operator’s  Conferences.  Dr.  Cochran  spoke  at  Purdue  University  and 
Dr.  Grayson  at  the  University  of  Massachusetts. 

Dr.  Grover  C.  Smart,  Jr.,  has  joined  the  staff  of  the  Holland  sub¬ 
station  as  Nematologist.  He  received  his  B.  A.  and  M.  A.  from  the  Uni¬ 
versity  of  Virginia  and  the  Ph.  D.  degree  from  the  University  of  Wis¬ 
consin. 

Several  staff  members  in  the  Agronomy  Department  at  Virginia 
Polytechnic  Institute  have  recently  received  grant-in-aids.  Dr.  G.  W. 
Thomas  will  study  “Electrolyte  Inhibition  by  Soil”  on  a  $19,000  NSF 
grant,  and  Drs.  R.  E.  Blaser  and  R.  E.  Schmidt  have  received  $2,000 
from  the  American  Potash  Institute  to  conduct  a  project  on  the  “Effect 
of  Potash  with  Rates  of  Nitrogen  and  Phosphorus  on  Tifgreen  Bermuda 
Grass  Turf  and  Winter  Growing  Ryegrass.”  Dr.  Blaser  has  also  received 
a  $500  grant  from  the  Golf  Course  Superintendents  Association  of 
American  for  an  undergraduate  scholarship  in  turf. 

New  members  of  the  Agronomy  Department,  Virginia  Polytechnic 
Institute  are  Mr.  H.  H.  Thomas,  Mr.  S.  J.  Ross,  Jr.,  and  Dr.  J.  Crews. 
Mr.  Thomas  and  Mr.  Ross  will  be  stationed  with  the  soil  survey  group 
in  Madison  and  Dr.  Crews  will  head  up  the  department’s  program  in 
flue  cured  tobacco  breeding. 

The  Virginia  Tech  Chapter  of  A.D.S.A.  (student  Dairy  Club)  was 
honored  as  the  outstanding  dairy  club  in  the  United  States  at  the  recent 
A.D.S.A.  meeting  in  Logan,  Utah. 
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